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Abstract. The recently developed method (Paper 1) enabling one to investigate the evolution of dynam¬ 
ical systems with an accuracy not dependent on time is developed further. The classes of dynamical 
systems which can be studied by that method are much extended, now including systems that are; (1) 
non-Hamiltonian, conservative; (2) Hamiltonian with time-dependent perturbation; (3) non-conservative 
(with dissipation). These systems cover various types of N-body gravitating systems of astrophysical and 
cosmological interest, such as the orbital evolution of planets, minor planets, artificial satellites due to 
tidal, non-tidal perturbations and thermal thrust, evolving close binary stellar systems, and the dynamics 
of accretion disks. 

PACS. 98.52.-b Stellar systems 


1 Introduction 

The typical procedure used in simulations of N-body gravitating systems [T] is the numerical iterative integration of 
the equations of motion, but this has a well-known difficulty: the inevitable accumulation of errors as the number of 
iterations increases and hence beginning from a distinct time further computations become meaningless. The difficulty 
is serious in the case of typical nonlinear and/or many dimensional systems covering, as now is understood, most of 
the physical systems. The systems of N gravitating bodies involving various astrophysical problems possess chaotic 
properties (e.g. [HinillllSlISMZ] ), which make their long term simulations rather complex. 

Recently a method of numerical investigation was developed [5] , providing a principal possibility for determination 
of the evolution of dynamical systems with an accuracy not dependent on time: an example of an evolution of each 
parameter of a system at given accuracy was exhibited in [S]. This method was here shown to be applicable for 
Hamiltonian systems with small periodic perturbation not dependent on time [S] 

H(J,z?,/?) = i7o(/)+/3Tfi(/,i?), (1) 

where Hq{I) is the integrable Hamiltonian, 1,^ are the action-angle variables and /3 is the small parameter of the 
non-integrable Hamiltonian Hi. 

Here we represent a further development of that method extending the class of dynamical systems which can be 
investigated, paying mainly attention to the systems having primary physical and astrophysical interest. We show that 
systems of rather general classes can be principally studied, including systems that are: 

1. non-Hamiltonian, conservative; 

2. non-conservative (with dissipation); 

3. Hamiltonian, with time-dependent perturbation. 

In our analysis we essentially use a result by Chernoff m concerning the properties of operators on a complete 
Riemannian manifold. 

Our initial aim is to investigate a rather general type first order differential equations 

= rix) ■ ( 2 ) 

The main idea of the method developed in [S] was the reduction of this equation to the following one 


i'll -I- Lu = 0 


(3) 
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where L is a differential operator. The solution of this equation by means of the computation of the resolvent using 
computer algebraic analytical codes gives the evolution of the function u{t) in time. Then the functions being 

determined by the function u(t) and describing the initial system, can be found (see [S])- 


2 The resolvent 

Let M be a d-dimensional complete Riemannian manifold with a metric g. Denote by F{M) the set of smooth functions 
determined on M. Let Fog{M) be the subset of F{M) which includes functions having a compact support. On Fq{M) 
one has the operation of inner product 


{u,v) = f u*{x)v{x){det gy^‘^d‘^x ; 

Jm 


(4) 


where u and v are elements of Fq{M). Denote by L'^{M) the completion of Fq{M) with respect to that product. 
Consider the operator 

L : F{M) F{M) 

in the form 


du 


Lu = -i'Y^ f°-{x)—a = -i^dau; 


where /“(a:) are smooth functions. 

One can show that the following relation for this operator is fulfilled 

{u,Lv)= f u*{x)Lv{x){det g)^^'^d‘’‘x 

Jm 

= i [ u* {x)v(x) {div f) {det g^^^d^^'x + {Lu,v), 

Jm 

divf = (detg)"^/^i9a((detg)^/^/“) = /“;a. 
divf = 0, 

(m, Lv) = {Lu, v ), 

i.e. L is a formal self-adjoint operator. 

In our case the velocity of the propagation c{x) introduced by Chernoff m is equal to 

c{x) = \\f{x)\\ = {gabrfy^'^ ■ 

According to Chernoff (Theorem 2.2 in [10]), if M is a complete Riemannian manifold and 


where u,v G Fo{M), 
When 
then 


/ CXD 

dr/c{r) = -Poo, 


(5) 


( 6 ) 


(7) 

( 8 ) 

(9) 


where c(r) = sup{c(a:);a: G S{y,r)}, S{y,r) is a ball of radius r with center y, then the operator L defined on space 
LJ{M) with domain Fq{M) is essentially self-adjoint. 

Consider now Eq.(3), to which, as mentioned above, Eq.(2) is reduced. 

The solution of that equation can be found by means of the Laplace transform |S], 

f d\ 

u{x,t) = - / exp{-iXt)Rx{l3)u{x,0)—■, (10) 

Jr 27ri 

F = {a + is , —oo < a < -Poo , s = sq > 0} (11) 

where R\{f3) is the resolvent, 

Rx{^) = {X-L)-\lmiX)>0. (12) 

Let us first consider the case when 
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r{x) = e+PC{x) 

and 

(13) 

div^ = 0 = div(, P = const ^ 1, 
then the resolvent is well approximated by formula [S] 

N 

(14) 

RxiP) = Y. p'^{Rx{0)B)'^Rx{0) + o{p ^), 

where 

Rx = {X-Lo)-\ 

Lo = -iCda , 

B = -iCda ■ 

In the case when for the dynamical system defined on M 

(15) 

divf ^ 0, 

we shall consider another one, having in local coordinates the form 

= r{x), 

y = -y divf. 

Here the manifold M is reduced to M x R with the metric 

(16) 

g + dy'^ . 

In the particular case that 

(17) 

M'" = Tor’^ X R^-’^ = Id-k} 

and 

r = C{I),a = l,...,k, 

!’> = 0,b = 1,... ,d - k , 

the resolvent Rx{P) can be found as in [S]. 

(18) 


3 Conclusions 


Let us summarize the dynamical systems covered by the above analysis. 
The systems represented by equations 


= nx) , 


include: 

1. Systems for which the Liouville theorem is fulfilled, 


divf = 0. 

This problem is reduced to finding out of the resolvent i?A(0), including the case when: 

a) the numbers of action-angle variables do not equal to each other. The calculations of Rx{0) should be considered 
separately at each given case. For example, R\{0) is easily found and therefore the method is trivially applicable when 

b) /“ is a periodic function; 

2. Systems with perturbed Hamiltonian 




(19) 


When H{I, d, t, /3) is a periodic function of time the resolvent is trivially represented via Fourier series (cf. [5]). In 
other cases the Rx{0) must be calculated separately; 
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3. Non-conservative systems, when 

divf ^ 0. 

This case is reduced to P.1, as shown above. 

The evolution of all these dynamical systems can be investigated by the method developed in [8] with an accuracy 
not dependent on the value of time. Evidently, now the classes of physical and astrophysical problems are far extended. 
Namely, besides traditional astrophysical nonlinear dynamical systems - N-body gravitating systems - now this will 
include an entire bunch of problems, e.g. the non-Hamiltonian motion of minor planets and space probes at tidal, non- 
tidal perturbation modes, and thermal thrust; these are particularly important, for example, for accurately testing 
of frame dragging, e.g. m- Other examples of non-conservative dissipative systems include e.g. the orbital evolution 
of close binary stellar systems due to evolving tidal effects, the magnetohydrodynamical evolution of non-stationary 
accretion disks near black holes, see e.g. UMMM- 
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